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Abstract 

An uncoupled multi-scale homogenization approach is used to estimate the effective 
thermal conductivities of plain weave C/C composites with a high degree of porosity. The 
geometrical complexity of the material system on individual scales is taken into account 
through the construction of a suitable representative volume element (RVE), a periodic unit 
cell, exploiting the information provided by the image analysis of a real composite system 
on every scale. Two different solution procedures are examined. The first one draws on the 
classical first order homogenization technique assuming steady state conditions and periodic 
distribution of the fluctuation part of the temperature field. The second approach is con- 
cerned with the solution of a transient flow problem. Although more complex, the latter 
approach allows for a detailed simulation of heat transfer in the porous system. Effective 
thermal conductivities of the laminate derived from both approaches through a consistent 
homogenization on individual scales are then compared with those obtained experimentally. 
A reasonably close agreement between individual results then promotes the use of the pro- 
posed multi-scale computational approach combined with the image analysis of real material 
systems. 

Keywords thermal properties, finite element analysis, micro-mechanics, carbon-carbon plain 
weave composite 
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(a) (b) 

Figure 1: Color images of a real composite system: (a) Scheme of multiscale structural model 
(from top - transverse and longitudinal view of fiber tow composite, composite unit cell, compos- 
ite lamina, composite plate), (b) Carbon tow microstructure showing major pores and transverse 
cracks. 



1 Introduction 

Carbon-carbon (C/C) plain weave fabric composites belong to an important class of high- 
temperature material systems. An exceptional thermal stability together with high resistance 
to thermal shocks or fracture due to rapid and strong changes in temperature have made these 
materials almost indispensable in a variety of engineering spheres including aeronautics, space 
and automobile industry. Applications include components in spacecraft protective shields, 
wing leading edges or parts of jet aircrafts turbine engines. 

While their appealing thermal properties such as low coefficients of thermal expansion and 
high thermal conductivities are known, their prediction from the properties supplied by the 
manufacturer for individual constituents is far from being trivial since these systems are gener- 
ally highly complicated. Apart from a characteristic three-dimensional (3D) structure of textile 
composites the geometrical complexity is further enhanced by the presence of various imperfec- 
tions in woven path developed during the manufacturing process. A route for incorporating at 
least the most severe imperfections in predictions of the mechanical properties of these systems 
has been outlined in [1] in the context of statistically equivalent periodic unit cell. Although 
properly accounting for three-dimensional effects the resulting representative volume element 
still suffers from the absence of the porous phase, which in real systems, Fig. 1, may exceed 
30% of the overall volume. As suggested in [2], neglecting the material porosity may severely 
overestimate the resulting thermal properties of textile composites. In this regard, the proce- 
dure introduced in the current paper extends the previous studies, as they either neglect the 
porous phase [3, 4] or analyze the porosity effects on the level of fibers and tows only [5]. 

The porosity of C/C composites directly arises as a result of the manufacturing process 
characterized by thermal decomposition and transformation of an initial polymeric precursor 
into the carbon matrix through several steps of carbonization, re-impregnation and final graphi- 
tization. As evident from Fig. 1 the major contribution to the porosity is due to crimp voids 
and delamination cracks, which are usually classified as inter-tow voids (Fig. 1(a)), as well as 
due to intra-tow voids represented by pores and transverse cracks developed within the fiber 
tow composite (Fig. 1(b)). For more details the reader is referred to [6] and references therein. 

While simplifying averaging schemes may provide rational micromechanics models for the 
prediction of effective thermal conductivities [7] at the level of the fiber tow composite in 
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Fig. 1(b), the complex mesoscopic structure of woven fabric plotted in Fig. 1(a) calls for con- 
siderably more accurate treatment of the actual geometry as already demonstrated in [1, 8]. 
In such a case, the image analysis combined with a reliable morphological description of the 
underlying composite structure then provide a general tool for the determination of what has 
been termed the statistically equivalent periodic unit cell (SEPUC) introduced by the authors 
in their previous works on random and imperfect composites [9, 1, 8]. Unlike classical averaging 
schemes, information on the local fields in representative volume elements formulated on the 
bases of SEPUCs is derived through a detailed numerical analysis which typically employs the 
finite element method (FEM). A special treatment of boundary conditions is then needed to 
establish a link with the first-order homogenization scheme which develops upon the assumption 
of homogeneous effective (macroscopic) fields [10]. 

Existence of macroscopically uniform fields (strains or stresses in the case of mechanical 
problem or temperature gradient and heat flux in the case of heat conduction problem) readily 
allows for splitting the local fields into macroscopic and fluctuation parts which in view of the 
solution of heat conduction problem reads 

9(x) = H x + 9*(x) or 6{x) = H iXi + 0*(x), (1) 

where H represents the macroscopically uniform temperature gradient vector and 0*(x) is the 
fluctuation part of the local temperature 9{x). Following, e.g. [11, 12] the solution of Eq. (1) then 
turns into the search for 6* in terms of the applied macroscopic uniform temperature gradient H 
or the macroscopic uniform heat flux Q. Consistency between the macroscopic (homogenized) 
quantities and volume averages of the corresponding local fields then requires either setting the 
boundary values of #* equal to zero or subjecting the fluctuation part of temperature field to 
periodic constrains. While both types of boundary conditions are equally applicable ensuring 
the macroscopic heat flux being equal to the volume averaged microscopic (local) heat flux, 
the latter conditions will be employed in Section 3 when developing the framework for heat 
conduction problems as they were shown to provide the best approximation for a fixed RVE 
size in a purely mechanical analysis, see e.g. [13, 14]. The need for periodic boundary conditions 
also arises when departing from asymptotic homogenization, see e.g. [15]. 

It is also interesting to point out that Eq. (1), when applied to either constrains on 
is consistent with so called afflne boundary conditions represented in our particular case by 
homogeneous temperature or flux applied on the outside boundary T of the RVE as 

9{x) H • x, Q n ( x ) = Q n on T, 

6{x) = HiXi, q n (x) = Qim, on T, (2) 

where n is the outer normal to T. Although the solution of a steady state heat conduction 
problem driven by prescribed macroscopic temperature gradient (or boundary temperatures 
consistent with Eq. (2)i), especially if the phase thermal conductivities are temperature inde- 
pendent, provides directly the effective conductivity matrix, it appears useful, particularly for 
more complex geometries as those in Fig. 1(a), to run the time dependent transient heat con- 
duction problem which allows for extracting considerably more information regarding thermal 
behavior of the composite material. Quite often, as will also be the case in this study, a com- 
mercial code, which does not allow for direct introduction of periodic boundary conditions, is 
used. In such a case, the boundary conditions given by Eq. (2) prove particularly useful as they 
naturally ensure that the volume average of microscopic (local) fields (temperature gradient or 
heat flux) are equal to their macroscopic (prescribed) counterparts H and Q as evident from 
Eqs. (9) and (11). 

To construct a certain periodic unit representing an entire laminate with all relevant geo- 
metrical details (distribution fibers within the fiber tow composite, waviness of fiber tow path, 
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porosity, etc.) might, however, prove rather impractical particularly from the computational 
point of view. Instead a so called uncoupled multi-scale approach [16], still at the forefront of 
material science interest, seems rather attractive allowing us to address the material complexity 
separately at different levels. Three particular levels of interest can be identified for the textile 
composite under consideration. Henceforth, they will be referred to as micro-scale (the level of 
individual fibers within a fiber tow, Fig. 1(b)), meso-scale (the level of a fiber tow composite, 
Fig. l(a)i_3 from top to the bottom) and macro-scale (the level of a laminated plate, Fig. l(a)4), 
respectively. To estimate a response of a such complex structure it appears reasonable to per- 
form a sequence of uncoupled analyses corresponding to individual scales. For this approach 
to be successfully utilized it is then viable to establish a link between individual scales. Here, 
the concept of scale separation plays a crucial role in the sense that the three analyses can be 
carried out independently such that output from one is used as an input to the other in terms of 
volume averages of the local fields while taking into account the boundary conditions mentioned 
in the above paragraphs. Such an approach is also adopted in the present study leading to a 
consistent search for the effective (macroscopic) thermal conductivities of a plain weave highly 
porous fabric composite. 

The paper is organized as follows. Following the introductory part our attention is paid in 
Section 2 to the formulation of various unit cells associated with individual scales. Theoretical 
formulation of the homogenization procedure for a steady state heat conduction problem is 
outlined in Section 3. Section 4 then illustrates the efficiency and reliability of the applied 
multi-scale analysis by comparing the numerical results with those derived experimentally. The 
essential findings are finally summarized in Section 5. 

2 Image analysis and construction of the geometrical model 

It has been demonstrated in our previous work, see e.g. [9, 1, 17, 18] that image analysis of real, 
rather then artificial, material systems plays an essential role in the derivation of a reliable and 
accurate computational model. This issue is revisited here for the case of woven fabric C/C 
laminate with particular relation to the adopted uncoupled multi-scale solution strategy. 



Table 1: Material parameters of individual phases [19, 20] 



Material 


Thermal conductivity 


Specific heat 


Mass density 




[Wm^K- 1 ] 


[Jkg^K" 1 ] 


[kgm -3 ] 


Carbon fibers 


(0.35, 0.35, 35) 


753 


1810 


Carbon matrix 


6.3 


1256 


1400 


Voids filled with air 


0.02 


1000 


1.3 



For illustration, let us now consider an eight-layer carbon-carbon composite laminate. In- 
dividual plies are made of plain weave carbon fabric Hexcel 1/1 embedded in a carbon matrix. 
Each filament (fiber tow) contains about 6000 carbon fibers T800H based on Polyacrylonitril 
precursor. Such fibers are known as having a relatively low orderliness of graphen planes on 
nano-scale. Nevertheless, unlike glassy carbon, they still posses a transverse isotropy with the 
value of longitudinal thermal conductivity considerably exceeding the one in the transverse di- 
rection, see Table 1. As already mentioned in the introductory part the carbon matrix forms 
as a result of several cycles of carbonization/densiflcation and final graphitization of the initial 
polymeric precursor (green composite) where the phenolic resin UMAFORM LE is used as the 
bonding agent (matrix). Note that phenolic resins belong to a family of non-graphitizing resins 
so that the final carbon matrix essentially complies, at least in terms of its structure, with the 
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Figure 2: Representative segment of eight-layer plain weave fabric laminate. 



original cross-linked polymeric precursor. Therefore, the resulting material symmetry is more 
or less isotropic with material parameters corresponding those of glassy carbon [21, 22]. 

A typical segment of the composite laminate appears in Fig. 2 showing characteristic porosity 
which may exceed at the structural level (macroscale) 30% [23, 24] and is often considered as 
an intrinsic property of this type of composite [21]. Several such micrographs were processed 
with the help of LUCIA G [25], Adobe Photoshop, Corel Draw and Matlab R12 softwares to 
acquire information regarding the basic structural units like an average thickness of carbon 
tows, size of voids, shape and essential dimensions of fiber tow cross-section, distribution of 
transverse and delamination cracks etc. which were subsequently exploited in the construction 
of representative unit cells on individual scales. 

2.1 Micro-scale 

Starting with the fiber tow composite as the basic structural element we recall Fig. 1 showing 
a typical shape of the fiber tow cross-section and significant amount of transverse cracks and 
voids resulting in a non- negligible porosity up to 15%. Unfortunately, a detailed analysis of 
the fiber tow cross-section would be computationally infeasible. As a suitable method of attack 
appears on the other hand formulation of a two-step homogenization problem. 

To proceed, let us consider a typical micrograph of the fiber matrix composite shown in 
Fig 3(b) taken as a random cut from the fiber tow cross-section evident in Fig. 3(a). Providing 
this section is sufficiently large to be statistically representative of a real microstructure it 
becomes possible to proceed in the footsteps of our previous work [9] and formulate a statistically 
equivalent periodic unit cell at the level of individual fibers. Keeping in mind on the other hand 
a relatively high volume fraction of fibers, approx. 55%, the results presented in [9] and remarks 
put forward in [26] allows us to conclude that the actual microstructure can be replaced with 
a simple periodic hexagonal unit cell plotted in Fig. 3(d). The resulting effective material 
parameters then serve as direct input for the analysis at the level of the fiber tow cross-section. 
Here, a new unit cell, again exploiting information acquired from image analysis, is introduced to 
properly account for the porous phase. The proposed periodic unit cell, which not only reflects 
the voids volume fraction but to some extent also their arrangement, is seen in Fig. 3(c). 

2.2 Meso-scale 

Having derived the effective material parameters for the fiber tow composite, the assumed three 
level homogenization procedure continues along the same lines on meso-scale. To that end, let 
us recall the representative section of the composite laminate in Fig. 2. A detailed inspection 
of this micrograph reveals three more or less periodically repeated geometries. For better view 
we refer to Figs. 4(a)-(c). 

Several such sections taken from various locations of the laminated plates were examined, 
again with the help of image analyzer LUCIA G, to obtain averages of various parameters 
including segment dimensions, fiber tow thickness, shape of the fiber tow cross-section also 



(a) 



(b) 




(c) 



(d) 



Figure 3: Homogenization on micro-scale: (a) fiber tow composite, (b) fiber-matrix composite, 
(c) voids-composite periodic unit cell, (d) fiber-matrix periodic unit cell 



position, size and location of large vacuoles. Approximately 100 measurements were carried out 
for each segment and subsequently utilized in the formulation of corresponding periodic unit 
cells displayed in Figs. 4(d)-(f). 

2.3 Macro-scale 

The final, clearly the most simple, step requires a construction of the homogeneous, thought not 
isotropic, laminated plate. The stacking sequence of individual periodic unit cells, which were 
introduced in the previous section, complies with that observed for the actual composite sample. 
Partially for the sake of simplicity, but also to be consistent with the analyses performed on 
lower scales, the periodic boundary conditions are considered even on the macro-scale. Providing 
we are interested only in the bulk response of the laminate thus ignoring detailed variation of 
the local fields, this assumption does not yield a significant error in the desired estimates of 
the macroscopic coefficients of thermal conductivities. Clear evidence is available in Section 4 
comparing the numerical and experimental results where the former ones are derived from the 
theoretical grounds presented in the next section. 

3 Theoretical formulation 

We now proceed to establish a framework for the determination of the effective thermal conduc- 
tivities regardless of the material system considered in the previous section providing the same 
boundary conditions are applied on each scale. An attentive reader will notice a close similarity 
with the derivation of the effective elastic material constants, see e.g. [11] for that matter. 
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(c) (f) 

Figure 4: Homogenization on meso-scale: (a)-(b) PUC1 representing carbon tow-carbon matrix 
composite, (c)-(d) PUC2 with vacuoles aligned with delamination cracks due to slip of textile 
plies, (e)-(f) PUC3 with extensive vacuoles representing the parts with textile reinforcement 
reduction due to bridging effect in the middle ply 



3.1 Governing equations 

Adhering to indicial notation with a — ^ and a i — representing the time and space 

at oxi 

derivatives, respectively, the simplest form of the balance equation reads 



(p'C-)0 + <& = O, 



(3) 



where p s is the mass density and Cp represents the specific heat of a given phase s. If referring 
to micro-scale, for example, the superscript s may first represent the fiber or the matrix phase 
and in the second step of homogenization, recall Section 2.1, it may be associated with the 
fiber-matrix composite and voids. The phase constitutive equations follow from the generalized 
version of Fourier's law and are provided by 



-y s h s 



hi 



in fl s 



(4) 



where hi = 9^, x is the thermal conductivity matrix [Wm _1 K _1 ] and if) = % _1 is the thermal 
resistivity matrix. Providing a heat flux, consistent with Eq. (2)2, is imposed over the entire 
boundary V we arrive, with the help of Eq. (4)i, at the following boundary condition 



on r. 



The weak form of Eqs. (3)- (5) is then given by 



/ 



8d{pC p )9 + Se ti XijO,j df2 + / 6&q n dT = 0. 



(5) 



(6) 



It is interesting to show that under the steady-state conditions (6 = 0). Eq. (6) is essentially 
equivalent to 

(MiXijOj) = SHiQi, (7) 

where (a) represents a volume average of a given quantity, i.e. (a) = y^y adtl. This becomes 
evident after introducing Eqs. (1) and (2)2 into the second term of Eq. (6) to get 



■^JsiHiXi + ^QjTijdr 



1 

W\ 

SHiQjSij = SHiQi. 



5HiQj J Xitij dr + 



66*Q j n j dr 



(8) 
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The integral J r 56*Qjrij dT disappears providing we set the fluctuation part of the temperature 
field #* equal to zero on T or impose the periodic boundary conditions (the same values of #* on 
opposite sides of a rectangular periodic unit cell). Either choice of variation of #* on T readily 
ensures that (hi) = Hi since 



(9) 



where integration by parts was used to transform the volume integral into the integral over the 
boundary T. As already mentioned in the introductory part, the periodic boundary conditions 
will be used in the actual numerical analysis. 

Note that Eq. (7) essentially resembles the Hill lemma in the context of pure thermo- 
mechanical problem [12]. Here it shows the consistency of the entropy change at the two 
associated scales (micro-meso, meso-macro) [27, 10]. 

The boundary conditions (2)2 further imply the volume average of the local heat flux be 
equal to the prescribed macroscopic heat flux Qi 

<«) = Qi- (10) 

This immediately follows under steady state conditions since q^i = so that qi = (qjXi)j. The 
volume average of the local quantity qi then gives, see also [10], 

1^1 (Qi) = / ftdfi = / (qjXijjdti 
Jn Jn 

= J qjfijXi dr = J q n Xi dT = J QjUjXidY 

= Qj [ XijdQ = \n\Qj8ij = \Q\Qi. (11) 

3.2 Effective conductivity and resistivity matrices 

To proceed, we limit our attention to steady state conditions and employ, in view of the forth- 
coming finite element formulation, the standard matrix notation; e.g. [28]. Then, under pure 
thermal loading consistent with the boundary conditions (2)i and taking into account the fact 
that variation of a prescribed quantity vanishes, we receive the following form of Eq. (7) 

m T bd {h}) = 0, (12) 



In the framework of the finite element method (EEM) the vector {h}, recall Eq. (1), is provided 

by 

{h} = {H} + [B]{0 d *}, (13) 

where [B] stores the derivatives of the shape functions and {9^} lists the nodal values of the 
fluctuation part of the temperature field. Substituting from Eq. (13) back into Eq. (12) yields 
the resulting system of equations 

[K] TO = {R}, (14) 

where 

[K] = / [B] T be] [B] cm, (15) 

{R} = - / [B] T [ X ]{H}df2 = -[S] T {H}. (16) 
Jn 
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Solving for {9^} from Eq. (14) gives 

R*} = -[K]- 1 [S] T {H} = -[G]{H}. (17) 

Next, introducing the vector {9^} back into Eq. (13) provides the volume average of the local 
heat flux in the form 

{Q} = <{q}> = - JL I [x] [ [I] - [B] [G] ] dO{H}. (18) 

I I «/ Q 

Finally, writing the macroscopic constitutive law in the form 

{Q} = -\ X ] hom {H}, (19) 
readily provides the homogenized effective conductivity matrix [x\ horn as 

[x]h0m = W\ I M[["]-[B][G]]dfi. (20) 

I I J £2 

In actual computations the coefficients of the effective conductivity matrix [x] h ° m are found 
as volume averages of the local fields from the solution of three successive steady state heat 
conduction problems. To that end, the periodic unit cell is loaded, in turn, by each of the 
two (2D) or three (3D) components of {H}, while the others vanish. The volume flux averages 
normalized with respect to {H} then furnish individual columns of [x\ hom ' . The required peri- 
odicity conditions (the same temperatures {9^} on opposite sides of the unit cell) are accounted 
for through multi-point constraints. In our particular case it suffice to assign the same code 
numbers to respective periodic pairs. 

The derivation of the effective resistivity matrix may proceed along the same lines providing 
the unit cell is loaded by the prescribed macroscopic uniform heat flux {Q}. In this particular 
case, the volume average of the local temperature gradient is not known a priori. Eq. (7) then 
yields two sets of governing equations for unknown nodal values of #* and volume average of 
the local temperature gradient ({h}) = {H} in the form 



[L] [S] 
[S] T [K] 



{hi i _ r -{Qi 
m j i {o} 



(21) 



where the matrices [S] and [K] were already introduced by Eqs. (15) and (16) and the matrix 
[L] is given by 

[L] = / [ X ] dO. (22) 
Combining the above equation together with the macroscopic constitutive law 



{H} = -M hom {Q}, (23) 

gives after some manipulations the searched homogenized effective resistivity matrix [ip\ hom in 
the form 

M W =[[L]-[S][K]- 1 [S] T ]" 1 . (24) 
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4 Results of the numerical analysis 



The present section summarizes the numerical results of the proposed three step uncoupled ho- 
mogenization scheme. On each scale the relevant periodic unit cell developed in Section 2 was 
discretized into finite elements. The geometrical complexity of individual computational mod- 
els together with the desired periodicity constraints led to extremely fine meshes. While this 
may seem irrelevant from the steady state conditions point of view, it suddenly becomes very 
important when running the transient heat conduction problem. This is also why the periodic 
boundary conditions were disregarded in the latter case. While the results are again presented 
separately for individual scales, it is worthwhile mentioning that the output in terms of the effec- 
tive conductivities derived at a lower scale served directly as the input for the up-scale analysis. 
Regardless of the type of analysis performed, steady state or transient conduction problem, 
only the boundary conditions of type (2)i were considered and the effective conductivities were 
found through the procedure described in Section 3.2. 



4.1 Micro-scale 




As suggested in Section 2.1 the effective conductivities of the fiber tow composite were 
derived in two steps. First, an intact carbon fiber-carbon matrix composite was considered. Two 
computational models displayed in Fig. 5 were examined. Both 2D and 3D analysis was carried 
out. A simple rule of mixture (RM) was used to estimate the effective thermal conductivities 
in the fiber direction (z-direction) for comparison with the general 3D analysis 

k% == kg ~\~ c kg , (25) 

where c^, c m , k(, k™ represent the fiber and matrix volume fractions and conductivities in the 
z-direction, respectively. Their values are listed in Table 1. The homogenized conductivities 
then appear in Table 2. 

Table 2: Effective thermal conductivities [Wm _1 K _1 ] - Intact fiber tow 



Geometry 


k x 


ky 


k z - 3DFEM 


k z - RM 


hexagonal 


2.17 


2.17 


21.93 


21.96 


square 


2.10 


2.10 


22.04 


21.96 



As expected, there is a minor difference in the results provided by both models. Never- 
theless, the effective values derived from the hexagonal array model were further employed in 
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the subsequent analysis step which allowed us to introduce the intra-tow voids into the intact 
but already homogeneous fiber tow composite. Again, two computational models evident from 
Fig. 6 were studied. 




Figure 6: Finite element mesh of fiber tow composite including voids: (a) real distribution of 
voids, (b) approximate distribution of voids 



Hereafter, we were concerned with the two-dimensional problem only. The corresponding 
results are available in Table 3. Apart from a steady state analysis a transient heat conduction 
problem goverened by Eq. (3) was addressed. Clearly, solving this problem then calls for the 
effective mass density and specific heat on the level of an intact fiber/matrix composite. For 
this purpose a logarithmic rule of mixture was adopted in the form 

lnp = c f \np f + c m lnp m , (26) 
\nC p = c f \nCl + c m \nC™. (27) 

The analysis was performed with the help of FEMlab commercial code using the Heat Transfer 
Mode [29] . To simulate a unidirectional flow the unit cell was heated on the one side while zero 
temperature was prescribed on the other side. The initial temperature was also assumed to be 
equal to zero. In addition, a zero flux boundary conditions were prescribed on the remaining 
sides to represent isolated surfaces. This clearly yields the mixed boundary conditions on T, 
which in turn may violate the conditions (9) and (11) if the periodic boundary conditions are 
disregarded, see the discussion in the following paragraphs. The conservation condition (7), 
however, holds even in this case. 

The thermal conductivities were then estimated from Eq. (19) after reaching the steady state 
conditions (time independent temperature profile). In such a case Eq. 3 reduces to q^i = thus 
naturally enforcing Eq. (11). If for example a unidirectional heat flow along the macroscopic x- 
axis is considered, then the component of the effective heat conductivity is, in view of Eq. (19), 
provided by 

hom _ (Ox) / 9 q\ 

Xx ~ ~WY (28) 

The Xy hom component of the effective conductivity matrix [x\ horn is derived analogously. If no 
action is taken this result corresponds to the assumption of t* = on T. 

The results in Table 3 suggest a minor difference between the steady state and transient 
heat analysis. Among others the lack of periodic boundary conditions, much coarser mesh [6] 
and truly approximate estimates of the homogenized mass density and specific heat appear as 
the most critical sources for the resulting difference. One may also expect this difference to 
grow when moving up the scales particularly if assuming a certain consistency of the three step 
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Table 3: Effective thermal conductivities [Wm 1 K 1 ] - Fiber tow with hexagonal arrangement 
of fibers including voids 



Analysis 


Steady state 


Transient 


Voids distribution 




ky 


k z - RM 




ky 


k z - RM 


real 
approximate 


0.97 
1.12 


1.70 
1.77 


19.01 
19.01 


1.20 


1.60 


19.10 



procedure in the sense that the results derived on a lower scale using one type of analysis are 
transferred to a higher scale where the same type of analysis is performed. 

4.2 Meso-scale 



... 











(a) 



■■ ' M * ■" mmn^gp* ^ 



(b) 



■ ■ ............. 



(c) 

Figure 7: Finite element meshes; (a) PUC1, (b) PUC2, (c) PUC3. 



At this level the carbon fiber tow is treated as a homogeneous phase with the effective 
material parameters derived from the numerical analysis on micro-scale. Since these are provide 
in the local (fiber) coordinate system, it is desirable, in order to account for the tow waviness, 
to transform the corresponding effective conductivity matrix into the global coordinate system 
as 

lx(x)}T = [T(a01 T be]]™ (x) [T(a01 , (29) 

where [T(x)] is the continuously varying transformation matrix. The matrix [x( x )fg W then 
enters Eq. (20) where applicable. The computational procedure is, nevertheless, identical to the 
one outlined in the previous section. With reference to Section 2.2 the three periodic unit cells in 
Fig. 7 were examined. Table 4 then provides the resulting homogenized thermal conductivities. 
To accept the notable difference in the results from the two types of distinct analyses we address 
the reader to the comments offered in the last paragraph of Section 4.1. 

Although more complicated, solving the transient heat conduction problem allows us to re- 
ceive further information regarding the response of a composite to thermal loading including a 
gradual evolution of the temperature profile and time lag to reach the steady state conditions. 
Fig. 8 shows a typical graphical output of the results provided by FEMlab. These plots, in 
particular, correspond to the onset of steady state conditions showing already more or less con- 
stant temperature gradient (linear variation of temperature) for all unit cells with corresponding 
distribution of heat fluxes. 
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Tirna 1 b; Surface: temperature (T> 

— i r 



Tima 1 a; Arrow: heat flu* proportional 




(c) 

Figure 8: Simulation of transient heat conduction problem - temperature profile and heat flux 
in the direction parallel to the composite plate: (a) PUC1, (b) PUC2, (c) PUC3 
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Table 4: Effective thermal conductivities [Wm 1 K x ] - Representative unit cells of textile 
composites 



Analysis 


Steady state 


Transient 


Cell 


fc-longitudinal 


^-transverse 


fc-longitudinal 


fc-transverse 


PUC 1 


9.46 


2.27 


9.10 


2.30 


PUC 2 


9.03 


1.47 


7.40 


1.80 


PUC 3 


7.29 


1.53 


6.60 


1.75 




(c) (d) 

Figure 9: Textile laminate with: (a)-(b) regular arrangement of plies, (c)-(d) irregular arrange- 
ment of plies. 



4.3 Macro-scale 

The macroscopic analysis of heat conduction problem represents the final step of the proposed 
multiscale approach. From the computational point of view it requires a relatively simple 
analysis of a three-layer laminate stacked from three periodic unit cells shown Fig. 7. However, 
the geometrical details of these unit cells are no longer relevant. Instead, they are treated 
as homogeneous blocks with the assigned homogenized (mesoscopic) properties. The resulting 
macroscopic properties are stored in Table 5. 

To validate these results an additional numerical analysis was performed. Here, all geomet- 
rical details of the laminated plate were taken into account by creating a multi-layered periodic 
unit cell exploiting the geometrical models developed for the mesoscopic analysis. The selected 
stacking sequence was motivated by images of the real laminate such as the one displayed in 
Fig. 2. Two such unit cells were constructed, see Fig. 9. While the first unit cell, Fig. 9(a), 
assumed an ideal stacking with regular arrangement of individual unit cells, the second one, 
Fig. 9(c), allowed for a mutual shift of individual plies to approximate the actual geometry 
more accurately. It is interesting to see, Table 5, that all types of analyses (laminate, PUC 
with regular and PUC with irregular arrangement of unit cells) essentially provide the same 
estimates of the macroscopic thermal conductivities. Such a conclusion clearly advocates the 
use of multiscale analysis at least in the present context of linear steady state or transient heat 
conduction problem. 

To judge the quality of any computational approach purely from numerical experiments 
seems, however, rather shallow. To enhanced credit of a numerical analysis it is therefore 
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Table 5: Effective thermal conductivities [Wm 1 K x ] - Laminated plate (The number in paren- 
theses indicates the difference between a numerical value and experimental data.) 



Analysis 


Steady state 


Transient 


Geometry 


fc-longitudinal 


/c-transverse 


fc-longitudinal 


/c-transverse 


Laminate 
Macro cell regular 
Macro cell irregular 

Measured 


8.47 (15.3%) 
8.65 (13.5%) 
8.60 (14.0%) 


1.66 (3.75%) 
1.68 (5.00%) 

1.67 (4.37%) 


8.30 (17.0%) 
10.00 


1.80 (12.5%) 
1.60 



desirable to compare the numerical results with those obtained experimentally. For this type 
of composite the results from experimental investigation of the thermophysical properties are 
available in [30]. For details on the pulse transient method used in this work together with 
computational models required to relate temperature to the generated heat pulse we refer the 
reader to the above paper. Some basic information can also be found in [6]. The experimentally 
determined macroscopic conductivities appear in Table 5. Taking into account the possible 
errors in the determination of phase material parameters (carbon fibers and carbon matrix) on 
the one hand and errors associated with the laboratory measurements on the other hand adds 
further confidence in the presented three-level uncoupled multiscale homogenization approach. 
In this context, the simplified two-dimensional approach (with the exception to the lowest 
scale) adopted for the solution of a generally three-dimensional problem, a natural component 
of laboratory measurements, should also be added to a list of sources of possible errors. It is 
worth noting even under these simplifications, the errors associated with the analysis compare 
well with the results of fully 3D analyses reported in [4]. Therefore, it is expected that a reliable 
three-dimensional geometrical model constructed on the mesoscopic level will further improve 
the predictive capabilities of the multiscale solution strategy. This is the topic of our current 
research. 

5 Conclusions 

Three levels of hierarchy are introduced in this contribution to derive the effective thermal 
conductivities of a plain weave textile laminate. Different resolution of microstructural details 
are considered on individual scales for the construction of an adequate representative unit cell. 
Such a unit cell arises as a result of elaborate evaluation of images of a real composite sample. 
The geometrical complexity of these types of composites are mainly responsible for a slow 
progress in the formulation of a generally three-dimensional unit cell. Possibility to enlighten 
this subject is given in [1] promoting the construction of such an RVE by matching statistical 
characteristics of both the real composite and RVE, at this step, however, in the absence of a 
porous phase. Owing to a significant contribution of this phase to the overall volume of the 
composite ruled the choice of a simplified two-dimensional analysis for the present study. 

Clearly, the theoretical formulation, here developed on the basis of the first-order homog- 
enization, is quite general and space invariant. Both transient and steady state conditions 
were examined in this study. When considering solely the temperature boundary conditions of 
type (2)i allows the solution to be accepted even if neglecting the periodic boundary conditions. 
Proper averaging relations then still provide the correct means for the scale transition. This has 
been exploited when running the transient heat conduction problem. The periodic constraints, 
on the other hand, were imposed for the estimates of effective thermal conductivities under 
the steady state conditions. In this case, too fine discretization needed for the straightforward 
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introduction of periodic boundary conditions was of minor concern. Both approaches have 
shown, however, their potential in the derivation of the desired effective (homogenized) thermal 
conductivities of highly complex plain weave textile composites through the application of fully 
uncoupled multiscale homogenization scheme. 
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